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In the recent years, remarkable attention is drawn to Stirling engine due to noticeable advantages, for 
instance a lot of resources such as biomass, fossil fuels and solar energy can be applied as heat source. 
Great numbers of studies are conducted on Stirling engines and non-ideal adiabatic method is one of 
them. In the present study, the efficiency and the power loss due to pressure drop into the heat exchang¬ 
ers are optimized for a Stirling system using non-ideal adiabatic analysis and the second-version Non- 
dominated Sorting Genetic Algorithm. The optimized answers are chosen from the results using three 
decision-making methods. The applied methods were compared at last and the best results were 
obtained for the technique for order preference by similarity to ideal solution decision making method. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Increasing the price of fossil fuel, environmental pollution, and 
noise pollution due to the engines using this fuel, has caused addi¬ 
tional motivation for research on other types of power generation. 
High thermal efficiency, minimal pollution, reliability, maximum 
utilization and clean combustion are the major features that we 
expect from the new engines. In the recent years, the Stirling en¬ 
gines have attracted a lot of attentions due to high capability 
[1,2 .A Stirling engine presents a reasonable theoretical efficiency 
which can be closer to the Carnot efficiency, comparing with other 
reciprocating thermal engines [3 . According to the Stirling en¬ 
gine’s ability to achieve the highest possible efficiency, this engine 
always is noticed for the researchers. With this regard, the first 
classical analysis of Stirling engine was done by Schmidt [4]. 
Schmidt assumed that the gas is isothermal in the expansion and 
the compression spaces. This assumption reduces the complexity 
of the compression and the expansion spaces with considering 
constant temperatures in the spaces. 

In the real engines, during the process, the working spaces tend 
to the adiabatic mode mostly. Based on this fact, Finkelstein carried 
out the first analysis of non-constant temperature of the Stirling 
engine in 1960. In the analysis presented by him, each part of 
the engine (cooler, heater, regenerator, expansion and compression 
spaces) was considered as a control volume and the conservation 
laws of mass and energy were analyzed using the equation of state 
[5]. Martini provided computer simulation of Stirling engine with 


* Corresponding author. Tel.: +98 9121947510; fax: +98 21 88617087. 
E-mail address: Akasa@ut.ac.ir (A. Kasaeian). 

0196-8904/$ - see front matter © 2014 Elsevier Ltd. All rights reserved. 

http://dx.doi.Org/10.1016/j.enconman.2014.01.022 


considering five working spaces for it [6]. In this model, pressure 
drop and heat loss were intended. Also, Urieli and Berchowitz used 
a computer code for solving differential equations using Runge- 
Kutta fourth-order method. They presented an adiabatic model 
for non-ideal mode in order to improve the numerical predictions. 
In this method, the concept of the non-ideal regenerator was real¬ 
ized [7]. 

Abbas et al. developed a non-ideal adiabatic method. In their 
work, the regenerator space was divided into two parts and this 
model included losses such as the shuttle loss and regenerator loss 
[8]. Strauss and Dobson developed other model using the quasi- 
adiabatic analysis in which the losses of the regenerator were con¬ 
sidered [9 . Granodos et al. developed the quasi-steady flow model 
to calculate the pressure drop through the heat exchanger chan¬ 
nels. In their model, the engine was divided into 19 parts which 
10 control volumes were dedicated to the regenerator [10]. Tlili 
et al. developed the second adiabatic model based on the 
quasi-steady flow. In their modeling approach, they considered 
the shuttle losses, the internal and external conduction losses in 
the regenerator, and the loss of energy that lead to pressure drop 
in the heat exchangers [11 . Parlak et al. carried out the thermal 
analysis on the Stirling engine with the Gamma structure. The 
quasi-steady flow analysis was performed to achieve more pre¬ 
cised results. The thermal efficiency reached to 25% using nitrogen 
as working fluid under pressure of 6.5 bar and the heat source 
temperature of 873 I< [12 . Considering the different methods of 
optimization in the energy issues, which have been carried out 
recently, it is worth to do multi-objective optimization of which, 
some are explained as followings. 

Different engineering problems such as skyline computation 
and vehicle routing issues have applied multi-objective 
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Nomenclature 



Qii 

heat transfer to heater 

M 

total mass of the working fluid 

Qk 

Q. 

heat transfer to cooler 
heat transfer to regenerator 

m 

mass of the working fluid in the different parts of the 
engine 

P 

mean effective pressure 

Re 

Reynolds number 

s 

stroke 

Pr 

Prandtl number 

fr 

engine’s frequency 

St 

Stanton number 

T h 

temperature of the working fluid of hot side 

m 

objectives function 

Th 

T c 

temperature of the heat source 
temperature of the heat sink 

X 

vector of decision variables 

w c 

work done by compression space 

Subscripts 

V Se 

expansion swept volume 

rh 

interference between regenerator and heater 

V Sc 

compression swept volume 

kr 

interference between cooler and regenerator 

V de 

expansion dead volume 

he 

interference between heater and expansion chamber 

V dc 

compression dead volume 

ck 

interference between compression chamber and cooler 

A 

heat transfer effective surface 

h 

heater 

A W g 

wetted area of the metal net of the regenerator 

k 

cooler 

/ h 

heater length 

r 

regenerator 

k 

cooler length 

c 

compression chamber 

lr 

regenerator length 

e 

expansion chamber 

V r 

heater volume 

j 

number of the inequality constraints 

V h 

V k 

cooler volume 
regenerator volume 

K 

number of the equality constraints 

f 

friction factor 

Greek 

letter 

dh 

d k 

heater hydraulic diameter 
cooler hydraulic diameter 

s 

p 

0 

effectiveness of the regenerator 

efficiency 

crank angle 

d r 

regenerator hydraulic diameter 

C v 

specific heat at constant volume 

y 

specific heat ratio 

Cp 

G 

k 

specific heat at constant pressure 
working fluid mass flow 
thermal conductivity 

p 

dynamic viscosity 


optimization [13-15 . The solution rout of the multi-objective 
optimization cases is a highly achieved target which requires the 
simultaneous satisfaction of a number of different and conflicting 
objectives. Evolutionary algorithms (EA) were applied in an at¬ 
tempt to stochastically solve problems of this generic class during 
the 18th century [16 . Inquirying a set of solutions of which each 
satisfies the objectives at an acceptable level without being over¬ 
come by any other solution, is a reasonable solving method to a 
multi-objective quandary [17]. A possibly uncountable collection 
of solutions so-called Pareto frontier is universally represented 
by the multi-objective optimization. The evaluated vectors of the 
Pareto frontier illustrate the best feasible trade-offs in the objective 
function space. In this case, multi-objective optimization of various 
thermodynamic and energy systems have been of high interest 
attention of researchers these days [18-25]. 

For considering all the above-mentioned issues, we have stud¬ 
ied two objective functions including the power loss due to pres¬ 
sure drop into the heat exchangers and the efficiency. In 
addition, the multi-objective optimization is conducted with four 
decision variables including the temperature of working fluid at 
hot side, frequency, the stroke and the mean effective pressure. 

Considering the parameters that have been neglected by most 
of the previous works on the performance optimization of Stirling 
engine such as the temperature of working fluid at hot side, fre¬ 
quency, the stroke and the mean effective pressure, we selected 
our approach for the current study. In the previous studies, the 
multi-optimization and adiabatic methods have not been linked 
together for Stirling engines and the adiabatic systems have only 
been utilized for modeling. In this work, by implementing the mul¬ 
ti-objective optimization algorithms, the power loss due to pres¬ 
sure drop into the heat exchangers and the efficiency are 
minimized and maximized, respectively. 


2. Theoretical model 

In the constant temperature of the Schmidt analysis, it is as¬ 
sumed that the compression and the expansion spaces during the 
engine performance remain under constant temperature condi¬ 
tions [26 . This assumption provides a contradiction that the hea¬ 
ter, the regenerator, and the cooler do not contribute in the any net 
heat transfer engine cycle. For this reason, they seem to be of none 
use. The real engine working spaces tend to the adiabatic mode 
during the processes. This assumption requires that the net heat 
transfer during a cycle is supplied by heat exchanger.An adiabatic 
model is shown in Fig. 1 for better understanding. This model is 
composed of five parts which are connected together sequentially. 
In this model, the working fluid in the cooler and the heater are at 
constant temperature conditions namely T k and T h , respectively. 
Also, the wire matrix regenerator and the embedded fluid are 
placed in a linear distribution. In this model, the assumption is that 
the working spaces are at adiabatic conditions. Also, it is assumed 
that there is no leakage of the working fluid and the pressure drop 
is not created in the fluid flow path. The performed work under 
these assumptions, is dependent on the changes in the volume of 
the compression and the expansion spaces and heat is transferred 
between the outside and the working fluid in the cooler and the 
heater which are Qk and Qh, respectively. 

The most common method to obtain the required differential 
equations and the heat transfer in the heater and the cooler, is 
writing the energy equation and the ideal gas equation for each 
control volume. The assumed control volume is shown in Fig. 2 

With getting differential of the mass balance for the assumed 
system [7], Eq. (1) is obtained. 

Dm c + Dm k + Dm, + Dm h + Dm e = 0 (1) 
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c k r h e 



Qt <& 




By taking the logarithm of both sides of the equation of state, 
and taking differential from it and assuming the constant temper¬ 
ature and volume of the heat exchanger, the following equation 
can be written as Eq. (2). 

(2) 

pm v y 

By inserting Eq. (2) in Eq. (1), we have: 

Dm c + Dm e +Dp (— + — + —) = 0 (3) 

\P P P J 

However, the equation of energy can be written for space com¬ 
pression as followings: 

DQ C + CpT ck gA ck = DW C + c v D(m c T c ) (4) 


Since the compression space is assumed to behave adiabati- 
cally, hence DQc = 0 and the done work is DW C = pDV c .According 
to the law of mass conservation, the working fluid accumulation 
rate (Dm c ) is equal to the input mass to the control volume (-gA ck ). 

c v T ck Dm c = pDV c + c v D(m c T c ) (5) 


Thus, the accumulation rates of the fluid in the expansion space 
and compression space are stated as: 


Dm r = 


Dm P = 


pDV c + % Dp 

(6) 

RT ck 

pDV e + 

RT ke 

(7) 


By inserting Eqs. (6) and (7) in Eq. (3), the desired equation 
based on the pressure differential is obtained. 


Dp 


—Yp( 


DVc DVe 


ck 


he 




+ 


Tie 


( 8 ) 


The changes in the volumes of the expansion and the compres¬ 
sion spaces are defined as follows: 


V e = 0.5 X Use X (1 + cos(0)) + U de 


(9) 


V c = 0.5 x V Se x (1 - cos(0)) + 0.5 x V Sc x (1 + cos(0 - 0)) 

+ V dc (10) 

Also, the temperatures of T C k and The are determined by consid¬ 
ering the fluid flow direction (gA). In fact, the direction of fluid flow 
can be assigned as followings: 


Dm = gA x - gA 0 = m x -m 0 


(ii) 


Using Eq. (11) for each of the sectors of Fig. 1, the following 
equation is obtained: 


m ck = -Dm c 
m kr = m ck - Dm k 
m rh = m kr - Dm r 
m he = m rh - Dm h 


( 12 ) 

(13) 

(14) 

(15) 


The total produced work by the engine can be writing as the 
sum of the work done in the compression and the expansion 
spaces: 

DW = p(DV c + DV e ) (16) 

By replacing energy equation, employing the aforementioned 
equations, and imposing some simplifications, we reach to more 
appropriate form of the energy equation. 


DQ. + (Cpli/Tii - c p T 0 m 0 ) = 


c p pDV + c v VDp 
R 


(17) 


Work is not done inside the heat exchanger because the relative 
volumes remain constant. The following equations for the heater, 
the cooler and the regenerator can be written: 


DQ k = CvV * Dp - c p (T ck m ck - T kr m kT ) 


(18) 
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DQr = v ^ P c p (T kr m kr T Ih m rh ) 

(19) 

DQh = * £ P c p( T rhm rh T he m he ) 

(20) 


2.2. The loss power due to pressure drops in heat exchangers 


In Stirling engine, the friction of the fluid is related to the flow of 
gas through the heat exchanger and the pressure drop and this 
problem leads to reducing the output power. The Pressure drops 
can be calculated using the friction factor coefficient related to 
the Reynolds number. The pressure drop through the regenerator 
is calculated as followings [11,28]: 



2x/x|ixV r xGxl r 

m r x d 2 



where /, V r , G, l r and d r are the friction factor correlation, the regen¬ 
erator volume, the working fluid mass flow, the regenerator length 
and the regenerator hydraulic diameter, respectively. The friction 
factor correlation has the following form for regenerator. 

/ = 54 + 1.43Re 0 78 (22) 

The pressure drop through the heater [28]: 

2 x/x (ix V h xGxl h 

ai H = - ~2 - ( Z: >) 

m.\i x dfo 


/ = 0.0791 Re 075 


The pressure drop through the cooler [28]: 



2 x / x jli x Vi ( x G x lk 
m k x dl 


(24) 

(25) 


The friction factor correlation of cooler is the same as what was 
applied for heater. 

The total pressure drop is calculated as: 
dP = dP R + dP H + dP K (26) 


By placing the total pressure drop in the following equation, the 
amount of lost power in the heat exchangers is obtained: 


power loss 




2.2. Energy loss due to external conductivity 


The energy loss due to the external conductivity is formed due 
to the non-ideal heat regenerating of the working fluid, within 
passing through the regenerator. In the non-ideal heat regenerator, 
the energy stored by the regenerator during the gas transmission 
from the expansion chamber to the compression chamber is not re¬ 
turned to the working fluid, completely [7 . In the case of ideal heat 
regenerator, the coefficient s is equal to one. Energy loss due to 
external conductivity would be written as: 

Q-rloss 0 — B) * (Q.rmax — Qrmin) (28) 


The ideal regenerator coefficient is obtained from the following 
equation [29]: 


NTU 
NTU + 1 



The effect of the non-ideal regenerator is defined as followings 
by using the number of the transfer units (NTU): 


NTU 


St x A 


Wg 


2A 



where A wg is the wetted area of the metal net of the regenerator in 
the exposure to the working fluid. Also, St is the Staunton number 
that can be considered as [30 : 


Cf 0.46 x Re -04 
“ Pr 

(31) 

where Prandtl number (Pr) is considered as 0.7 for the utilized 
working fluid. The actual thermal energy in the cooler and the hea¬ 
ter spaces is calculated by the following relation: 

Qh Q-h "T Qrloss 

(32) 

0-K Qk — Q.rloss 

(33) 


Finally, the actual work of the Stirling engine is calculated by 
obtaining of the loss work due to the pressure drop in the heat 
exchangers. The actual work and the efficiency of Stirling engine 
are defined as below: 

actual work = (W - dW loss ) (34) 

actual heat = QH 
actual work 

n =- 

actual heat 

A method for simultaneous solving of these equations is using 
their initial values inside them. Thus, the solution process starts in 
a full cycle having the initial values of the all variables and integrat¬ 
ing them. This solution is carried out using the fourth order Rung- 
Kutta. The implemented engine for this research was the GPU-3 
whose characterizations are demonstrated in fable 1 . The GPU 3 en¬ 
gine is of the Beta type Stirling engine which has been described in 
details in the literature [31,32 . For the purpose of optimization, the 
theory of multi-objective optimization is described in next chapter. 

3. Multi-objective optimization 

In the multi-objective problems, due to a conflict between the 
goals, there is no single optimal solution. Thus, the method should 
search several optimal solutions that are called the Pareto set 
[33,34]. In the Pareto optimal series, by moving from one point 
(solution) to another, an objective function gets better and another 
objective function gets worse; because both solutions of the series 
are non-dominated relative to each other. 

The general form of multi-objective problems is as follows [34]: 

' Minimize/Maximize/ m (x) m = 1,2,..., M 

Subject to gj(x) ss 0 j = 1,2,... ,J 
< h k (x) = 0 fo = 1,2,... ,/C (i ’ 

,X[ L) Xj <x! U) i = 1,2,...,n 


(35) 

(36) 


Table 1 

The GPU3 engine parameters. 


Parameter 

Value 

Heat source T H (K) 

1200 

Heat sink T c (K) 

288 

Engine’s frequency (Hz) 

41.7 

Mean pressure P (MPa) 

4.13 

Regenerator length (mm) 

22.6 

Wire diameter d w (pm) 

40 

Regenerator diameter(mm) 

22.6 

Porosity 

0.697 

Thermal conductivity (Wm -1 K -1 ) 

16.6 

Total mass of gas (kg) 

0.001135 

Displacer diameter (mm) 

69 

Number of regenerators per cylinder 

8 

Number of gauzes of the matrix 

308 

Working fluid 

Helium 
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An x solution is a vector of n decision variables: 

X = [X 1 X 2 ...X n ] r (38) 

Each decision variable is limited between the amounts of the 
lower limit (Xj L) ) and the upper limit (Xj U) ). These limits make 
the space of the decision variable. Also,/is a vector of m objective 
functions. 

/(x) = [fl(x)/2(x).../n(x)] T (39) 

Of which gj (x) and h k (x) are constraint function that J and I< are 
number of the inequality constraints and number of the equality 
constraints, respectively. Although the single-objective optimiza¬ 
tion problems may have a unique optimal solution, the multi¬ 
objective optimization problems usually have a possibly uncount¬ 
able set of solutions on a Pareto front. Each solution associated 
with a point on the Pareto front is a vector whose components rep¬ 
resent trade-offs in the decision space or Pareto solution space. The 
multi-objective optimization problems evaluation function, /: 
Q -► A maps decision variables (x = Xi, • • • , x n ) to vectors (y = ai, 

• • •, a k ). This situation is represented in Fig. 3 for the case n = 2, 
m = 0, and k = 3. This mapping may or may not be onto some region 
of objective function space, dependent upon the functions and con¬ 
straints composing the particular multi-objective optimization. 
[35] 

3.1. NSGA-II 

NSGA-II is the second version of the “Non-dominated Sorting 
Genetic Algorithm” for solving non-convex and non-smooth single 
and multi-objective optimization problems. The concepts of 



Decision Variable Space 


Fig. 3. A schematic of the decision variable 


dominance and crowding distance are important subjects in the 
NSGA-II which were widely described by Deb et al. [37,38] 

Fig. 4 shows the performance of NSGA-II multi-objective opti¬ 
mization model. In this method, the children population (Qt) is 
made up by the parent population (P t ). Then two populations com¬ 
bine with each other and (P t ) population of size 2N is composed. 
Then, non-dominating sorting is used to classify the entire of the 
(R t ) population. Following the creating of different non-domination 
ranks (Fi, F 2 , • • •, Fj), the next population is filled in order of the pri¬ 
ority of these ranks. The next population (P t+ i) starts to fill the best 
non-domination ranks, then the second non-domination ranks and 
so on, until the P t+1 is filled. About the last rank, the solutions, 
which have more crowding distance, are the priority to fill P t+ 1 . Fi¬ 
nally, the population of children (Qt+i) from the parent population 
(P t+ i) is created using a select tournament function and the opera¬ 
tors of crossover and mutation. 

3.2. Objective functions, decision variables and constraints 

A mathematical model of the optimization process has three 
general parts (the objective functions, decision variables and con¬ 
straints problem) that are described here in. Two objective func¬ 
tions for this study are the power loss due to pressure drop 
inside the heat exchangers and the efficiency of the Stirling engine. 
In this paper, four decision variables have been considered as 
followings: 

/-: engine’s frequency, 
p: mean effective pressure, 

S: stroke, 

T h : the temperature of working fluid at hot side. 



space and objective function space [36]. 


f. 


q. 


Non-dominating sorting 


Crowding distance sorting 




Fig. 4. The structure of NSGA-II multi-objective optimization model [36]. 
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TOPSIS 

solution 


"1-T 

+ 


Non-ideal 


\ 
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390 400 

Power losses (W) 


410 


420 


430 



Fig. 5. (a and b) Pareto optimal frontier in the objectives’ space. 


The objective functions with respect to the following con¬ 
straints have been solved: 


' 40 < fr ^ 50 Hz 
4.1 ^ p ^ 6.89 MPa 
0.03 < S s: 0.07 m 
„ 900 < Th < 1100K 



For choosing the optimum point of the engine, a decision should 
be taken and there are different ways to make decisions. In this 
work, three methods including TOPSIS, LINMAP, and Fuzzy have 
been used which are explained briefly as followings. In the LINMAP 
method, the optimal point is a point on the Pareto Front that has 
the shortest distance to the ideal point. The ideal point is a point 
that has best situation. In the TOPSIS method, the optimal point 
is a point which has the lowest value of C L . 



dr 

d i+ + dr 



where d [+ and d r are the distance of each point from the positive 
ideal point and negative ideal point, respectively. In the Fuzzy 
method, the Fuzzy Bellman-Zadeh is used which is explained in fur¬ 
ther details in the literature [39-42 . 


4. Results and discussion 

The power loss due to pressure drop into the heat exchangers 
and the efficiency of Stirling engine are minimized and maximized, 


respectively. These objectives were maximized simultaneously 
using a class of multi-objective evolutionary algorithms called 
NSGA-II. The optimization process was performed using objective 
functions expressed by Eqs. (27) and (36) along with the con¬ 
straints which are expressed in Eq. (40). Fig. 5a illustrates the Par¬ 
eto frontier in the proposed objectives’ space obtained in the 
optimization scenario. Three final solutions were selected by the 
Fuzzy Bellman-Zadeh, LINMAP and TOPSIS decision-making meth¬ 
ods which are indicated in this figure. According to Fig. 5a, the ob¬ 
tained points by LINMAP and TOPSIS are approached towards each 
other. Also, it was shown that the optimal value of the power loss 
varied from 365 W to 425 W and the optimal value of t] was be¬ 
tween 0.475 and 0.52. Based on the obtained curve fit depicted 
in Fig. 6b, Eq. (42) was developed. 

f] = a x exp(b x powerloss) + c x exp(d x powerloss) (42) 

The coefficients (with 95% confidence bounds) and goodness of 
fit are summarized in fable 2. 


Table 2 

The coefficients and goodness of fit. 

Coefficient 

a = 19.33 b = -0.00631 c = -55.22 d =-0.009958 

Goodness of fit 

SSE: R-square: Adjusted R-square: RMSE: 

0.0001045 0.9709 0.9611 0.003407 















P (bar) f r (HZ) 
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Fig. 6. The distribution of f T for the optimal points on Pareto front. 
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Fig. 7. The distribution of P for the optimal points on Pareto front. 
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Fig. 8. The distribution of T h for the optimal points on Pareto front. 
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Fig. 9. The distribution of S p for the optimal points on Pareto front. 



Table 3 

Comparison of the final optimal solutions of this paper (multi objective-multi variable) with corresponding results. 


Method 

T h (K) 

S (m) 

P (bar) 

fr (Hz) 

Power loss (W) 

ri(%) 

Initial theoretical status 

939.30 

4.6 

41.3 

41.7 

441.759 

48.73 

TOPSIS 

1069.1 

3.03 

49.78 

40.73 

375.1 

49.65 

LINMAP 

903.9 

3.00 

41.32 

40.03 

382.4 

50.42 

Fuzzy 

1059.2 

3.01 

47.88 

40.53 

384.8 

50.75 

Urieli and Berchowitz [7] 

- 

4.6 

41.3 

41.7 

- 

62.06 

Timoumi dynamic model without losses [43] 

- 

4.6 

41.3 

41.7 

- 

54.96 

Timoumi dynamic model with loss dissipation (Ml) [44] 

- 

4.6 

41.3 

41.7 

- 

53.10 

Timoumi dynamic best model [44] 

- 

4.6 

41.3 

41.7 

- 

52.50 

Experiment [44] 

- 

4.6 

41.3 

41.7 

- 

35.00 


Table 4 

Error analysis based on the mean absolute percent error (MAPE) method. 

Decision-making method 

Fuzzy 


TOPSIS 


LINMAP 


Objectives 

Power loss 

4 

Power loss 

4 

Power loss 

4 

Average error (%) 

17.16 

1.92 

2.44 

0.06 

4.90 

2.10 


Figs. 6-9 demonstrate the scattered distribution of the design 
variables versus population. These figures can provide a good vi¬ 
sion on the variation of the decision variables from the Pareto fron¬ 
tier. Distribution of variables for the optimum points on the Pareto 
front for f r over population are exhibited in Fig. 6 and the range of f r 
varied from 40 to 40.8. It can be seen from Fig. 7 that the magni¬ 
tude of P for optimal Pareto front varied from 41 to 50 bar. As 
shown in Fig. 8, the variation of T h for the optimal Pareto front var¬ 
ied from 900 I< to 1080 K. Finally, the variables distributions for the 
optimum points on Pareto front for S p over the population are 
exhibited in Fig. 9. The upper and lower ranges of S p are positioned 
from 3.05 to 3.00. 

The optimal results for the decision parameters and objective 
functions using LINMAP, TOPSIS and Bellman-Zadeh decision-mak¬ 
ing methods are summarized in Table 2. The comparison of the 
optimal results obtained in this paper using three decision making 
methods with the corresponding results obtained in Refs. [7,43,44] 
is demonstrated in Table 3. Also, the initial status in where the 
parameters had not been optimized yet, are reported in the table. 
In this state, the code was written and the model was obtained, 
without considering optimization. With multi-objective optimiza¬ 
tion, we could decrease the power loss caused by the pressure drop 
in the system. 


According to Table 3, the characteristics of this work are com¬ 
pared to the other works which have not considered the important 
losses. This reason caused that they reported higher efficiency. 
Also, if the efficiency is optimized as single-objective, the efficiency 
will be higher than the reported efficiency while this work has 
been carried out as two-objective optimization. 

The mean error analyses for the objective functions are catego¬ 
rized in Table 4. The analysis is done for all mentioned decision¬ 
making methods. 

5. Conclusions 

In this study, adiabatic analysis has been applied to determine 
the power loss and the efficiency of a Stirling engine. The power 
loss and the efficiency of the Stirling engine were considered 
simultaneously for multi-objective optimization and the engine’s 
frequency (f T ), mean effective pressure (p), stroke (S p ) and the tem¬ 
perature of the working fluid of hot side (T h ) were considered as 
the design parameters. Multi-objective evolutionary algorithm is 
considered based on NSGA-II algorithm and the Pareto optimal 
frontier in the objectives space was obtained. A final optimal solu¬ 
tion was selected from the solutions of the Pareto frontier using 
three decision making methods including LINMAP, TOPSIS and 
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Fuzzy methods. It was shown that the multi-objective multi-vari¬ 
able approach leads to more desirable design of the Stirling engine 
if it is compared with the single-objective single-variable ap¬ 
proaches which have been performed in the previous researches. 
Also, with comparing the results, the TOPSIS decision making 
method leaded to a better conclusion in aspect of power loss. 
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